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1— H Abstract. Among various definitions of quantum correlations, quantum discord 

has attracted considerable attention. The connection between the quantum 

I I discord and the entanglement of formation is described by Koashi- Winter relation. 

„^i We investigate this relation from the viewpoint of the quantum channel that is 

C^H isomorphic to the given state. It is shown that in the case of two-qubit states the 

channel, on the one hand, determines the form of the quantum steering ellipsoid 

of the given state, and on the other hand, is closely related to the concurrence of 

jrt the complement state of the given state. We also point out that, for two-qubit 

f—^ rank-two state, von Neumann measurement is the optimal choice to achieve the 

«H quantum discord. However, for some two-qubit states with the rank larger than 

i i two, the three-element POVM measurement is optimal. 
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1. Introduction 

Quantum correlations between parties of a bipartite quantum system have been the 
significant and interesting subject in quantum information processing as well as in 
quantum theory. There are various measures to quantify the quantum correlations, 
among which quantum discord has attracted considerable attention. 

Quantum discord is originally defined as the the difference between two classically 
identical (but quantumly distinct) formulas that measure the amount of mutual 
information of a pair of quantum systems [1] . An alternative way to define quantum 
discord is by writing it as the difference between total correlation and classical 
correlation [2], The total correlation is usually quantified by the mutual information 

l = S(p A ) + S( P B )-S(p AB ), 

where p is the state of a bipartite quantum system composed of particle A and 
particle B, the p A ^ = Tcb(a)(p AB ) is the reduced state of A (B), and S(p) is 
the von Neumann entropy of a quantum state p. As for the classical correlation, 
we adopt the definition given by Henderson and Vedral [2], Suppose that a positive 
operator valued measure (POVM) measurement is performed on particle A. The 
set of POVM elements is denoted by M = {M k } with M k > and J2k M k = 1- 
The operation element M k will produce particular outcome, labeled by k, with the 
probability pk — Tr{p AB (M k <8> !)]• When the outcome k occurs, the postmeasurement 
state of particle B is p B = Tr A [p AB (M k <g> l)]/p k . That is, the POVM measurement 
on A will "realize" a postmeasurement ensemble {p k ,p k } for particle B. In other 
words, the reduced state p B is decomposed into J2 k VkP B due to the measurement M. 
performed on A. The amount of information acquired about particle B is then given 
by S(p B ) — ^2 k PkS(p B ), which depends on measurement Ai. This dependence can 
be removed by doing maximization over all POVM measurements, which gives rise to 
the definition of classical correlation: 
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Denote the average entropy ^2 k PkS(p k ) as S , and the minimal value as S n 
Quantum discord is defined as 

D A->B =I C A^B = s{pA) + s A ^ B S(p A B) (2) 

The key step, and also the major difficulty, is to work out the minimal average entropy 
(MAE). So we usually focus on the MAE 5 mln . We call the POVM measurement 
on A optimal if such a measurement generates the MAE for B. The optimal POVM 
may not be unique. 

Similarly one can define classical correlation C Ai ~ B and quantum discord D A ^ B , 
corresponding to the case that the measurement is performed on particle B to gain the 
information about particle A. Generally, the quantum discord and classical correlation 
defined from A to B may not be equal to that defined from B to A, meaning that 
they are asymmetric. In this paper, we only consider the direction from A to B. For 
the sake of clarity, we will omit the superscript "A — > B" in the expressions of average 
entropy, classical correlation and quantum discord. 
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The concept of quantum discord is proposed to describe the quantum correlations 
which are not limited to entanglement. There is a fundamental difference between 
entanglement and discord for mixed states, although they are equivalent for pure 
states. A typical example is that there exist non-entangled (or separable) states 
with nonvanishing discord [1]. Quantum discord has received much attention 
in studies involving fuzzy measurement [3], mixed-state quantum computation 
speedups [4, 5] , broadcasting of quantum state [6] , complete positivity of dynamics 
[7, 8] , complementarity and monogamy relationship between classical and quantum 
correlations [9-11], dynamics of discord [12,13], operational interpretations of quantum 
discord in terms of state merging [14, 15], and the relation between discord and 
entanglement [16-19]. 

Let us focus our attention on the Koashi- Winter (K-W) relation presented in [11]. 
There the authors show that there exists a monogamy relation between the classical 
correlation of a bipartite state p AB and the entanglement of formation (EoF) of 
the complement p BC . By "complement", we means that both p AB and p BC are 
the reduced states of a tripartite pure state p ABC , namely, p AB = Trc p ABC and 
pBC = Tta p ABC_ Then K _ w relation tells us that C(p AB ) +E(p BC ) = S(p B ), where 
C(p AB ) is the classical correlation (from A to B) of p AB and E(p BC ) is the EoF of 
p BC . By noting the definition (1), we can rewrite K-W relation as S m [ n = E(p BC ). 
Obviously, K-W relation establishes the relationship between the EoF and the classical 
correlation (or minimal average entropy, or quantum discord) in an extended Hilbcrt 
space. It is remarkable that the particle B is the pivotal point of the relationship. 

In this paper, the pivotal role played by particle B is investigated from the 
viewpoint of quantum channel. It is well known that there is an isomorphism between 
quantum states and quantum channel [20-24]. Restricted to a two-qubit state p AB , 
we can say that, associated with p , there is a channel A acting upon the qubit B 
of the Bell state \$> AB ) = ^ 5 (|0 A )|0 B > + |1 A )|1 S )). We further show that that the 
significance of the quantum channel A is twofold. On the one hand, A determines the 
quantum steering ellipsoid [25], which has been used to discuss the quantum discord 
in [26]. As such, we present an alternative approach to the quantum steering ellipsoid. 
On the other hand, A is closely related to the concurrence of the complement state 
p BC . In fact, since p BC is a rank-two state, its concurrence can be evaluated exactly 
by means of the methods provided in [27, 28] . 

Subsequently, we point out that K-W relation implies the following two facts, 
(i) For any 2-qubit rank-two state, the quantum discord or classical correlation can 
be evaluated exactly by using K-W relation. In this case, we need only to perform 
POVM measurement with two (rather than three or more) operation elements, i.e., 
von Neumann measurement, on qubit A to achieve the optimal two-state ensemble for 
p B . The two states in the optimal ensemble have the identical von Neumann entropy. 
In addition, we present the geometric description for the optimal ensemble in terms of 
the quantum steering ellipsoid, (ii) For some 2-qubit states with rank larger than two, 
the von Neumann measurement may not be optimal, while the three-element POVM 
measurement will give a larger classical correlation or a smaller quantum discord. This 
fact has been implied by K-W relation, together with the numerical results obtained 
in [28]. 

This paper is arranged in the following manner. In Section 2, we consider the 
isomorphism between quantum states and quantum channels, and show that for two- 
qubit states the quantum channel determines the quantum steering ellipsoid. It is 
found that the ellipsoid is invariant under any filtering operation performed on a 
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specific qubit. In Section 3, a given two-qubit state is purified to a tripartite pure state, 
and the concrete process is presented to evaluate the concurrence of the complement 
state. We find that the quantum channel obtained in Section 2 is most relevant to 
the evaluation of the concurrence. Section 4 is a discussion about K-W relation. 
There we indicate that the von Neumann measurement may not be optimal to give 
the quantum discord or classical correlation for general two-qubit state. Section 5 is 
devoted to exact expressions as well as the geometric picture for the quantum discord 
and classical correlation of any 2-qubit rank-2 states. In addition, the conjecture 
proposed in [26] is proved. Two-qubit states with rank larger than two are discussed 
in Section 6. We show that some known results about the EoF can lead further insight 
into the quantum discord and classical correlation. 

2. Quantum channel and quantum steering ellipsoid 

It has been shown that quantum steering ellipsoid is a useful tool to discuss quantum 
discord and classical correlation [26]. In this section, we present an alternative 
approach to establish this geometric picture, that is, we will show that the quantum 
steering ellipsoid eventually boils down to the notation of quantum channel. 

2.1. Quantum, channels, one-qubit case 

A quantum channel $ is a trace-preserving, completely positive map. We refer the 
reader to textbooks (e.g. [29] and [30]) for a detailed description of quantum channel. 
Let p be the density operator on n-dimcnsional Hilbert space. The action of <I> on a 
state p can be expressed in the operator-sum representation: 

p^ $[p] = Y, K iP K l 

i 

with Ki the Kraus operators and J2i K\Ki = 1. We arc able to represent the channel 
$ in the form of superoperator, i.e., a n 2 x n 2 matrix: $ = J^-Ki ® K* . With the 
state p represented as a "long" vector in n 2 -dimensional Hilbert space, denoted by 
\p)), we can write $[p] as &\p))- In the following, we consider the simple one-qubit 
case. 

Let X be any operator on the two-dimensional Hilbert space, namely, 



Define the 4-vector (in column form) X as X = {xqq,Xqx,X\q, £ii) T ! where the 
superscript T means the transpose. The channel $ is now expressed as a 4 x 4 
matrix, and the action on X is given byX^ $X. On the other hand, the operator 
X can also be written as X = ^ 5Zu=o x \*.°~\n with o~o the 2x2 identity matrix 1 



2 



and Oi for i — 1, 2, 3 the Pauli matrices. That is, X is expressed by another 4-vector 
x = (xq,xi,X2,x^) . Define a unitary matrix T, which is from [31] 

/ 1 1 \ 

110 

i -i [6) 

\ 1 -1 ) 

We see that x = ^/2TX. Then in terms of the 4-vector x, the action of $ is 

x -> x' = T$T f a; = L^x, (4) 



T 1 

V2 
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where Lq, = T$T^. 

By unitary operations, L§, can be expressed in the canonical form [32]: 

Mil)- (5) 

where k = (ki, K2, k^) t is a column vector in three-dimensional real space, and S is a 
3x3 diagonal matrix, 5 = diag(£i,£2i£3)- For one-qubit state p—\ Eu^o^c' 7 *" the 
4-vector x = (x ,xi,X2,x 3 ) T is transformed to y — L$x. In this sense, we call L$ 
the Bloch representation of the channel $. 

Notice that x is distributed in the Bloch ball (including the surface), namely, 
x\ + x\ + x 2 ^ 1, or equivalently, x T r\x ^ with 77 = diag(l, — 1, — 1, —1). It follows 
that y is constrained by 

y T L-JnL^y > 0. (6) 

Here it is assumed that L$ is invertible. It is not difficult to see that (6) indicates an 
ellipsoidal region. Then the geometric meaning of (5) is clear: i$ transforms the unit 
sphere S 2 in M 3 (i.e., x 2 + y 2 + z 2 = 1) to an ellipsoidal surface, given by 

{x- Kl ) 2 (y-K 2 ) 2 (z-k 3 ) 2 _ 

tl ' c2 ' p2 L - 

SI ?2 ?3 

2.2. Quantum channels, two-qubit case 

Now we consider the effects of quantum channels on the two-qubit states. A 2-qubit 
state p AB can be expressed in Hilbert-Schmidt space as p AB — \ ^ V _ Q R^o-^ <g> a v 
with i? M „ — Tr(p AB a ll ® <r v ). Define i?-matrix as R = (i? M „). The matrix R can 
be also expressed as R = 2T(p AB ) R T T , where (p AB ) R is the reshuffling of p AB [31]. 
Given a bipartite matrix 

X = Yl VWJ'lvWfl 
the reshuffling of X is defined as 

x R = £ x iti/;jtjl \ij){i'j'\. 

Now suppose that qubits A and B are acted upon by local channels <& A and $ B 
respectively. The state p AB is mapped to p' — ($ A <g) <& B )[p AB ]. Expressing <& A 
and $ s as 4x 4 matrices, we have the following formula, which is essentially from [30]. 

{p AB )R _^ {p ,AB )R = q> A {p AB ) R^B ) T_ (?) 

Expressing it in the Bloch representation, we have 

R — > R 1 = L A R(L B ) T , (8) 

where R' is the i?-matrix of p' AB , L A = T$ A Tt and L B = T$ B Tt. 
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2.3. Quantum steering ellipsoid 

It is well-known that there is an isomorphism between quantum states and quantum 
operations, in particular, quantum channels [20-24]. In this subsection, we use this 
isomorphism to derive the equation of quantum steering ellipsoid from the quantum 
channel. 

Given a 2-qubit state p AB , we assume that the reduced state p A = Tr^ p AB is 
invertible. Define the state p AB as 

P AB = ^=P AB -^. (9) 

Then p A = Tib p AB = |-"-2- The isomorphism between state and channel is expressed 
as 

p AB = (id®A){P + ], (10) 

where id denotes identity operator and P + = |$ + )($ + | with |$ + ) = |(|00) + 1 1 1) ) a 
Bell state. 

According to formula (7), we have (p AB ) R = P R A T . By noting that P R = gJ.4 
with 1 4 the 4x4 identity matrix, we see that A T = 2(p AB ) R . In the Bloch 
representation, the channel is represented as 

L A = TAT f = 2T[(p AS ) fi ] T T t 

= 2[T*T^T(p AB ) R T T } T = R T R T + , (11) 

where R is the i?-matrix of p AB and R+ is the i?-matrix of P+, that is, R+ = 
diag(l,l,-l,l)=T*Tt. 

According to (6), the output states of the channel A are constrained in an 
ellipsoidal region, which is given by 

y T [L-JnLl 1 ] y = y T [R^nR^] y ^ 0. (12) 

Here we assume that La (or R) is invertible. Define the matrix $7 = R T r]R. It is the 
matrix f£ -1 that determines the ellipsoid, which we denote by l£. 

Define a similar matrix SI = R T nR with R being the i?-matrix of p AB , i.e., 
R = 2T(p AB ) R T T . An ellipsoid <B can be defined from il~ x . The equation is given by 

y T [R' l r,R' T ] y = 0. (13) 

We will show that £ is identical to <£. 

In fact, the ellipsoid € (or (5) is invariant under any local filtering operation 
performed on qubit A. To see this, let us consider how the matrix Q is transformed 
under the local filtering operation. The process of any local filtering on A is given by 

p AB -^(F®l)p AB {Fi®l), (14) 

where F is a 2 x 2 nonsingular matrix with F'F ^ 1. The right-hand side of (14) is 
an unnormalized state. By defining $j? = F <g> F* , we rewrite (14) as 

(p AB ) R — ► $ F (p AB ) R . (15) 

Note that $i? is not a channel, since it is completely positive but not a trace-preserving 
map. In the Bloch representation, (15) is equivalent to R — > LpR with Lp — T<&i?TT. 
By noting that Lp^nLp = \ det(_F)| 2 7y, we have 

R T nR — ► R T L T F i 1 L F R = \ det(F)\ 2 R T r 1 R. 
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That is, ft — > | det(F)\ 2 fl for any local filtering operator F. It follows that the ellipsoid 
(£ is invariant under the operation F (g> 1. 

It is remarkable that (13) is just the equation of the quantum steering ellipsoid 
which we introduce in [26] to discuss the quantum discord. Then we arrive at the 
following conclusion. 

Proposition 1 Given a two-qubit state p AB , assume that the matrix R and the 
reduced density matrix p are all invertible. Let p = (2p ) _1 ' 2 p (2p )~ 1 ' 2 . 
Then the state p is isomorphic to the channel A, which is expressed in Bloch 
representation as L\ = R T R^_, where R and R+ are, respectively, the R-matrix of 
the p AB and of the Bell state -4=(|00) + |11)). The channel La, in turn, determines 

the quantum steering ellipsoid (£, which is given by y T L~^ r\L~^ y — 0, or equivalently, 
y T Rr 1 r\R^ T y — 0. Moreover, the ellipsoid 2; is invariant under any local filtering 
operation on qubit A. 

We add a remark here. In the above discussion, the channel A (or La), which is 
isomorphic to the state p , is applied on qubit B. If we are concerned with qubit 
A, the quantum steering ellipsoid will be determined by the matrix R~ T rjR~ 1 , just as 
we presented in [26]. 

3. Quantum channel and concurrence of complement state 

In this section, we will show that the channel L\ is closely related to the concurrence 
of the complement state of p AB . According to [11], the state p BC is complement to 
p AB when there exists a tripartite pure state p ABC such that Tr a(p ABC ) = p BC and 
Tr c (p ABC )=p AB . 

Assume that the rank of p AB is four. The eigendecomposition of p AB is 

3 

i=0 

We need a four-level quantum system, denoted by qudit C, to purify p AB . The 
purification of p AB is 

i*> = E^i^i ?c >' ( 16 ) 

i=0 

where \i c ) for i = 0, 1, 2, 3 constitute the basis of the Hilbert space of qudit C. Then 
pBC _ Tr ^ |v|/)(\l/| i s a 2 ® 4 rank-2 state and is the complement to p AB . 

Note that (16) can be regarded as the Schmidt decomposition with respect to the 
bipartite partition AB : C. We can also perform Schmidt decomposition with respect 
to the partition A : BC. Assume that the eigendecomposition of p BC is 

P BC = ao|<po)v?o| +ai|Vi)(¥'i|- (17) 

It follows that 



|^) = v / "ol Q; o>!S)|¥'o)+v / ^ri a i>®l^i> ; (18) 

where \oti) for i = 0, 1 are the eigenvectors of p . In fact en are the eigenvalues of p 
associated with |atj). 

To proceed to evaluate the concurrence of p BC , we need the following two 
definitions, which are from [27]. 
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Definition 1 For any nxn matrix A, the second symmetric function 82(A) is defined 
as 

S2(A) = ^[(Tr A) 2 -Tr A 2 }. (19) 

Let H d be the space of complex Hermitian matrices of size d x d and let H+ 
be the cone of positive semidefinite matrices. The concurrence of a bipartite state is 
given by 

Definition 2 Given a d\®d 2 bipartite state p, denote by pi and p 2 the reduced state of 
the first and the second subsystem, respectively, i.e., Pi(2) — Ti"2(i) p. The concurrence 
Con(p) for the state p is defined on H + x 2 as the convex roof of the function 2y // S 2 (pi) 
or2y/S 2 {p2)- 

Now let us consider the concurrence of the 2 ® 4 bipartite state p B . The first 
step is to write out the reduced state p B . Let us consider a more general case: an 
arbitrary Hermitian operator X BC on the two-dimension Hilbert space T-l 2 spanned 
by the two eigenvectors of p BC , i.e., FL 2 = span{\ipo), \ifii)}- The reduced operator 
X B = Trc X BC . Define the following four operators on T-L 2 . 

?o = IvoX^oi + l^iX^il, <^i = l^oX^il + IViX^ol, (20a) 

? 2 = -i\ip )(ipo\ +i\ i Pi)( l Po\, ?3 = \<Po)(ipo\ - IviXvil- ( 20& ) 

These four operators are analogous to the identity matrix cto and the Pauli matrices 
Ui with i = 1,2, 3. We can call them generalized Pauli operators. Then we can express 
X BC z& 

X BC = \J2^ (21) 

with x^ real numbers. 

To express X B explicitly we have to calculate Trc \fii)(fv\ for p, v = 0, • • • , 3. 
For this purpose, we resort to a maximally entangled state of the bipartite system 
A : BC, that is, 

l*> = ^(|0 A )®|^> + |l A )®|^i)). (22) 

The state |\P) can be obtained from |^) given by (18) by applying local filtering on 
qubit A, followed by an unitary operation: 

|*) = [07F)®1 4 ]|*), 
where F — 2 — |cto)(c*o| + -7k=\ a i) { a i\i an d the unitary U transforms |ao) and |aj.) 
to \0 A ) and \1 A ), respectively. Now we can see that 

for i, j = 0, 1. It follows that 

Trc\v i ){p j \=2{i A \[Tr c \*)^\j A )=2{i A \p AB \j A ), 

with p AB — Trc I^X^I- Considering the generalized Pauli operators given by (20a) 
(206), we have 

Tr c?0 = 2Tr A [(a <g> l 2 )p AB ], Tr c ?1 = 2Tr j4 [(a 1 <g> l 2 )p AB ], 
Tr c ^2 = -2Tr A [(a 2 <g> l 2 )p AB ], Tr c c 3 = 2Tr A [(a 3 <g> l 2 )p AB }. 
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Now we can express X B as 

X B =Tr c X BC = l 2 J2 x ^c^ 



= Tr A { [(ocoo-o +x 1 a 1 - x 2 o- 2 + x 3 cr 3 )<E>l} p AB } . (23) 

Straightforward calculation shows that 

X B = ±(x T R + R) ■ a = \{x T Ll) • <r, (24) 

where er = (cto,<ti,<72, 0-3), and in the last equation we have used (11). 

In order to obtain the concurrence of p B , we use the following theorem, which 
is essentially from [27]. 

Theorem 1 Let X be a d\® d 2 bipartite matrix of rank not exceeding 2. Further, let 
H2 C <C dld2 be a linear complex subspace of dimension 2 such that the range of X is 
contained in R.2- Denote by H dld2 the space of all Hermitian matrices on C dld2 . Let 
H2 be the subspace of all matrices in H dld " 2 whose range is contained inR.2- Define the 
two quadratic forms Qi : A^ 2[(TrA) 2 - Tr(Tr 2 A) 2 ] and Q 2 : A^ (TrA) 2 - TrA 2 
onH d ^ 2 . 

Then the generalized eigenvalues of the pencil Q\\h 2 ~ wQ2\h 2 are oil real. Denote 
them by W\, W2, w 3 and W4 in decreasing order. Then the concurrence of X is given 
by 

Con(A) = y/Qi(X)-w 2 Q 2 (X). (25) 

In the case we are considering, Qi(X) is given by 

Q^X 80 ) = 2[(Tr X BC ) 2 - Tr(Tr c X BC ) 2 ] 

= 2(TrA B ) 2 -2Tr(A s ) 2 

= 2x — (x L A ) 

= x T L A nL A x, 

where in the third equation we have used (24). Then the matrix form of Q\ is 

Q 1 = L T A nL K . (26) 

As for Q 2 (X), we have Q2(X BC ) — ^x t tjx. It follows that the matrix form of Q2 is 
given by Q 2 = \r}. 

According to Theorem 1, we have to solve the following equation to find w 2 . 

det [Ll v L A - |t?) = 0. (27) 

It is equivalent to 

det [kriR T - ^r?) = 0. (28) 

Or in other words, w is the eigenvalue of 2RrjR T -q. 

Obviously, the channel L A plays the central role in the evaluation of Con(p BC ). 
It should be mentioned that one can employ the approach presented in [28] to get the 
same results. 
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4. Purification and Koashi- Winter relation 

Given a bipartite state p AB , by denning its complement state p BC , K-W relation [11] 
establishes the connection between the classical correlation of the former (i.e., p AB ) 
and the EoF of the latter (i.e., p BG ). This relation gives new insight into both of the 
optimal ensemble with respect to the classical correlation or the quantum discord and 
that with respect to the EoF. We outline briefly the K-W relation and make some 
remarks as follows. 

We restrict our attention to the case of two-qubit states. Given a p AB with rank 
r (r = 2, 3, 4), we add a r-dimensional ancilla C. The purification of p AB is then given 

by 

r-l 

l*)=E^®l^ B )K C ) ; (29) 

t=0 

where the probabilities qi and the pure states \ipf B ) constitute an ensemble for p AB , 
and {|i )}»=o,---,r-i is the basis of the Hilbcrt space of particle C. Note that \ip AB ) 
may not be the eigenstates of p AB and may not be orthogonal to one another, and 
that the rank-two reduced state p BC = Tr^ | v l')(5'| is the complement state of p AB . 

Suppose that an n-element POVM measurement A4 = {-^fe} with 
k = 0, • • • , n — 1 is performed on qubit A. We need only consider the case that all M& 
are of rank one [33,34], that is, M k is proportional to the one-dimensional projector. 
From the viewpoint of the whole system in the pure state | <I>) , the measurement A4 
gives rise to an ensemble for p BC , denoted by S BG — {p k , \f BC )}, where 

p k = (*|(Af fc ®l 2 ®lr)|*), 

Wk C )(fk C \ = Tr.4 [(M fc ® l a ® l r ) |*)<*|]/j» fcl 

with l r being r x r identity matrix. From the viewpoint of the state p , the 
measurement Ai on A gives rise to an ensemble for p B , that is, S B = {pk,P B }i 
where 

p k = Tr[(M k ®l 2 )p AB ], 



p B =Tr A [(M k ®l 2 )p AB ]/p, 



k- 

It is easy to see that the ensemble S B can be induced from S BC ', namely, p B = 
Tr c (\<pf c )(f B c \). In figure 1, we illustrate the relations just mentioned. 

The EoF of each pure \ip B c ) is given by E(ip BC ) — S(p B ). The average EoF over 
the ensemble <§ is given by 

n— 1 n—X 

E BC = Y,P k E^ B k C ) = Y,P^{p B k )- 
fe=0 k=o 

It follows that the average entropy of the ensemble /§ B is equal to the average EoF: 
S = E BC . (30) 

Note that (30) depends on the measurement on A. The dependence is removed 
by by performing minimization over all POVM measurements. It then follows that 
Smm = -^min- Koashi and Winter pointed out that the -E m ; n obtained in this way is 
exactly the EoF of p BC . Then the K-W relation is expressed as 

Smin = E{p BC ). (31) 
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It follows that 

C = S( P B )-E(p BC ), 



„AB\ 



D = S{p A ) + E(p^) - S{p AB ) 
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(32) 
(33) 



P i= 



M 



8 B ={P„ Pi) 



Purification 



Tracing out C 



|¥>0F| 



M A 



Tracing out C 



=K \<Pt c KK'\} 



Figure 1. The 2-qubit state p is purified into a tripartite pure state |\P). For 
p , performing a POVM measurement M on qubit A induces an ensemble for 
qubit B, namely <S B . For \^), the M B induces the ensemble for the composite 
system BC, namely S BC . Tracing out particle C from each member in S BC will 
give rise to S B . 

Let's make some remarks as follow. 

(i) Conventionally the EoF is regarded as coming from mathematical 
consideration. Now we see from K-W relation that EoF has relevance to the 
measurements. 

(ii) If rank(p" 4s ) = 2, then p BC is a two-qubit rank-2 state. It follows from [35] 
that E(p BC ) is achieved on the 2-state ensemble S BC such that E(p BC ) = E(ip BC ) = 
E(ipf c ). This means that we need only perform the two-element POVM measurement, 
namely, von Neumann measurement, on qubit A in order to obtain Smin and thereby 
the discord D. Moreover, the optimal measurement will induce an equi-entropy 
decomposition of p B , namely, S{p B ) = S(p B ). 

(iii) In the case of vank{p AB ) — 3 or 4, although the state p BC is of rank two, 
there is no definite answer as to whether a two-state ensemble will yield the EoF of 
p BC . To put it in terms of the discord D (or the MAE 5 m ; n ), there is no definite 
answer as to whether the optimal measurement performed on qubit A is the von 
Neumann measurement. Therefore, strictly speaking, the MAE obtained by means of 
von Neumann measurement only provides the upper bound of the EoF of p . 



5. Two-qubit rank-two states 

According to the above discussion, we see that the discord of any two-qubit rank-2 
state can be evaluated exactly. In this section, we provide the analytic results and the 
geometric description. 
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5.1. Analytic results 

The eigendecomposition of p AB is 

/ S = Ao|^o)(^o|+Ai|^i)(M (34) 

From the conclusion in [36], we know that the eigenstates \ipo) and \tpi) can be 
simultaneously transformed to the following forms by local unitary operations. 

|Vo>=ao|0)|0} + 6 |»?)|l), |Vi)=ai|l)|0)+6 1 |7; ± )|l), (35) 

where \ak\ 2 + \bk\ 2 — 1 for k = 0, 1, and |?7) is orthogonal to \r) ). Let 

\r,) = c\Q)+d\l), \ V x ) = -d*\0)+c*\l), (36) 

with |c| 2 + |d| 2 = 1. 

Attaching a qubit C to the two qubits A and B, we write the purification of p AB 

as 

I*) = V%|Vt>> <8> |0) + v/ATlV'i) <8> |1). 

The reduced state p BC is a two-qubit rank-2 state, and the concurrence can be 
obtained easily. The square of the concurrence is given by 

[Con(^)] 2 



— 2AoAi 
— 2AoAi 
The EoF of p BC is then 



\a b lC * - ai b c\ 2 + 2|d| 2 (|a | 2 |& 1 | 2 + K| 2 N 2 ) 
(ao&ic* — aiboc) 2 — 4aoaibobi\d\ 2 



£( ,-) = »(i±^HK>E), (37) 

where the function H(x) is defined as 

H(x) = -s;log 2 a;- (l-.x)log 2 (l - x), (38) 

for x € [0, 1]. From (32) and (33) we obtain the analytic expressions for the classical 
correlation C and quantum discord D. 

For later reference, we need the expression for [Con(p BC )] 2 in the case that the 
parameters a^, bk, c and d are assumed to be real numbers. In this case, we have the 
following results. 

(i) If c 2 (a 6i — ai&o) 2 ^ 4a ai6o6id 2 , then 

[Gon {i) (p BC )} 2 = 4A A 1 d 2 (a & 1 + M ) 2 . (39) 

(ii) If c 2 (aobi — aibo) 2 < 4aoaibobid 2 , then 

[Con (ii) (p BC )] 2 = 4A A 1 (a 6 1 - ai 6 ) 2 . (40) 

From (37) we have the corresponding E^{p BC ) and E^(p BC ). From (32) and (33), 
we have Cm, CVjn and -Dm, Dt^y 



Quantum discord of two-qubit rank-two states 13 

5.2. Geometric picture for entangled states 

As stated in Section 4, for any two-qubit rank-2 state, the optimal measurement on 
qubit A will realize an equi-entropy decomposition of the reduced state of qubit B. We 
will illustrate this result in the geometric picture presented in [26]. In this subsection 
we discuss the case of entangled states. 

For the entangled state p AB , the i?-matrix of p AB is invertible and thus the 
quantum steering ellipsoid l£ exists. Its equation has been given in Section 2.3, where 
the consideration is based on the notation of quantum channel. As a comparison, let 
us spend a moment to give a brief outline of the original discussion presented in [25] . 

Suppose that a POVM measurement is performed on qubit A. The 
postmeasurement state of qubit B, corresponding to the operation element Mk, is 
p B . In Hilbert-Schmidt space, we express M k , p A as 
3 1 3 

M k = Y^ *fc,/*o>. Pk = 2 Yl ^>M°M> 

H=0 fJ.=Q 



x k ,n = r Tr(M fc er p ), y k>fl = Tr(pf cr p 



where 

1 

2 
We define two 4-component vectors (in column form) x k and y k as 

Xk = (Xkfi, Xk,l, Xk,2, X k ^) T , 

Vk = (yk,o, Vk,i, Uk.2, Vk,z) T ■ 

Note that y k o = 1 for all k and the 3-component vector -ffk = (yk,i, Vk 2i Vk,3) T is just 
the Bloch vector of p B . Direct calculation shows the following equations. 

3 

Pk Vk = R T x k , Pk = ^ R ^° Xfe ^' ( 41 ) 

(41) provides the relationship between the measurement on A and the induced 
component in postmeasurement ensemble for the reduced state of B. Recall that 
R is invertible. So the vector t/k is in a one-to-one correspondence to the vector Xk- 
It follows that Xk = PkR^ T Vk- The non-negativity of Mj. leads us to 

vl[BT l riRr T ]v k >Q, (42) 

which is just the ellipsoid equation given by (12) or (13). 

The steering ellipsoid yields concrete geometric picture when we minimize the 
average entropy S. Denote by f B = (rf , r B , r B ) T the Bloch vector of p B . The 
state p B can be represented by a point B with the coordinate given by {rf , r B , r B ). 
It can be seen that point B is in the interior of £. If qubit A is measured by a 
n-element POVM measurement, the n Bloch vectors f B (k — 0, • • • ,n — 1) of the 
postmeasurement states p B constitute a convex polytope V, which is contained in the 
ellipsoid l£, that is, B G V C £. The average entropy is given by 

with r^ = |r^| and the function H given by (38). 

By noting that both V and £ are convex and the entropy function is concave, we 
see that the minimal value of S must be attained on the surface of (£. The surface 
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of l£ represents all postmeasurement states of qubit B when qubit A is measured by 
rank-one POVMs. Denoting by d€ the surface of £, we can write 

S min = min5> fc S(pf ) = min ^.//(li-lf), (43) 

where the minimization is taken over all {pk, f B } with r B £ <9(£ and ^^Pk^ 3 = T?S - 

Now we proceed to give the geometric picture. With p AB given by (34) - (36), 
we let all parameters a^, b^, c and d take real numbers in order to simplify our 
calculations. In this case the classical correlation and quantum discord have been 
obtained in Section 5.1. We present below the geometric description of the optimal 
decomposition of p B . 

Firstly let's define the following functions for later convenience. 

/i = A ao& + Aiai6i, f 2 = A a &o - Aiai&i, 

/3 = \AoAi(ai&o + a &i), / 4 = VAoAiOify) - a 6i), 

/s = v/A Ai(a ai + boh), f 6 = v/AoAi(a ai - b a b{). 

From (42) the equation of steering ellipsoid l£ is given by 

(c 2 /3 + d 2 fl )y\ + (fi + d 2 f 2 )y\ + 2cf 3 f 6 y iy3 

+ (c 2 f! + d 2 fl)y 2 2 - IchhVx - 2f 5 f 6 y 3 + (f 2 - d 2 f 2 ) = (44) 

In order to describe <£ more clearly, we apply a rotation about y 2 axis: 

Vi^yiCOsn-y 3 smi], y 3 -> y x sin 77 + y 3 cos 77, (45) 

with the rotation angle 77 given by 

sm?7 = =, cos?? = — . (46) 

Denote by <£' the ellipsoid after the rotation. The equation of €' is 

d'fiyl + ^fi + d 2 .!^ 2 



+ (/l + c 2 /! + d 2 A 2 ) 



„ Wg+gj 

^/f(/;:-/-:-^/f+d 2 /i a ) 



f2 + *ft + <pft (47) 

It can be seen that the 2/3 axis is one of the symmetry axes of £' and the center of (£' 
is located on the point 



0, 0, 



hVfi + c 2 fj 
fl + c 2 fl + d 2 f! 



Before rotation, the Bloch vector of the local state p B is 

f B = 2efiii + [\o(a 2 - bl) + X^al - b\)]e 3 , 

where ik denotes the unit vector along y^ axis. After the rotation, the Bloch vector 
f B is transformed to 

cos 77 sin 77 

r' u = I 1 1 r ls 

- sin -q cos 77 
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Figure 2. Schematic illustration of steering ellipsoid £' for p . Point B' 
corresponds to vector r' B , which is the Bloch vector of p B after rotation given 
by (45). A plane containing point B' and parallel to J/11/2 plane intersects the 
surface of l£' in a ellipse. Depending on the parameters of p , the optimal 
ensemble corresponds to either of the following two decompositions of r ,B : the 
convex combination of rjj and f^, or the convex combination of fp and rn. 



with the the angle 77 given by (46). Obviously, r' B is still in the 2/12/3 plane. 

In hgure 2, we illustrate the post-rotation ellipsoid (£' and the post-rotation Bloch 
vector r' B . Generally, the ellipsoid €' may not possess the rotational symmetry about 
any principle axis, and the Bloch vector r may not lie along the 2/3 axis. Nonetheless, 
there exist two forms of equi-entropy decomposition of f' B , which we describe as 
follows. 

Suppose a plane is parallel to j/ij/2 plane and passes through point B' . This plane 
intersects the surface of (£' , and the the intersection is an ellipse. Two chords of the 
ellipse, MB' N and PB'Q, pass through point B' . The former is parallel to y\ axis 
and the latter parallel to 2/2 axis. We see that OM = ON and OP — OQ. Then the 
equi-entropy decompositions are given by 

r' B = Pm r M +PNr N , (48) 

r' B =p P r P + p Q f Q , (49) 

where Tm denotes the vector from point O to M, and so forth. It follows that the 
average entropy corresponding to the decomposition (48) is given by 



>Smjv 
Similarly, (49) gives 



PmH 



l + r M 



PnH 



(- 



r N 



H 



l + r M 



S 



PQ 



H 



l + r P 



After tedious calculations, we get the length of OM and OP: 



' M 






- 4A Ai(i 2 (a &i + ai6i) 2 , 
4A Ai(a 6i - ai&o) 2 - 



(50) 
(51) 
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Recalling (37), (39) and (40), we have 

S M N = E {i) (p BC ), S PQ = E {ii) (p BC ). (52) 

Hence we have shown that (48) or (49), together with (50) and (51), provide the 
geometric description of the optimal postmeasurement ensemble for p B . 

5.3. Geometric picture for separable states 

If a two-qubit state p AB is separable, its steering ellipsoid may not exist. Nonetheless, 
the geometric description of the MAE is very clear, which we have presented in [26]. 
There we conjectured that equi-entropy decomposition is the necessary condition 
which must be satisfied in order that the average entropy reaches the minimum. Now 
the discussion in Section 4 shows that this conjecture is true. A complete description 
is given below. 

Any two-qubit separable state with rank two can be expressed, up to local unitary 
operations, as 

p AB = #)<0| ® |0)(0| + (1 - q)Wl^\ ® |</>)(</>|, (53) 

where q g (0,1), \ip) = cosa|0) + sina|l), and |0) = cos/3|0) + sin/3|l) with 
a, 8 G (0, J]. Here neither a or 8 takes the value of zero, otherwise the state p AB is 
a trivial product state. The purification of p is 



I*} =^|0)® |0) <g> |0) + ^T~q |V) ® 10) ® |1). (54) 

The concurrence of p BC is easily calculated, that is, 2^/q(l — p) cos a sin 8. From K-W 
relation (31), we have 



__ /l + a/1 - 4g(l - o) cos 2 a sin 8 \ 

S min = H(-^ - J (55) 

Then quantum discord D follows easily. 

Now we present the geometric description of the MAE. Suppose that we perform 
von Neumann measurement on qubit A. The measurement operators are given by 

M± = - ± xia x ± x 2 o- y ± x 3 a z , (56) 

with x\ + x\ + x\ = 1/4. Inserting the 4-component vectors x± = 

(l/2, ±xi, ±X2, ±£3) into (41) will give p± and y±. The concrete forms of p± 
and y± are a bit lengthy and we do not write them out explicitly. It can be shown 
that the Bloch vectors y± satisfy 

V±.3 +y±,itan/3 = 1. (57) 

Then the postmeasurement states of qubit B is constrained on the line L, given by 
(57). The reduced state p B is also on this line (see figure 3), as the Bloch vector of 
p B is given by 

= f(l- 9 )sin2/3, 0, q+{l-q)cos2f3) , 



f B 



which satisfies r B + r B tan/3 = 1. 

Now we consider the equi-entropy decomposition of p B , which means that 

r B = P+ y+ + P- y- with \y + 1 = \y_ \ . 
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Figure 3. The postmeasurement states of qubit B, as well as the local state p B 
itself, is constrained on the line L, the equation of which is given by (57). The 
optimal ensemble for qubit B corresponds to the case in which the length of OE 
is equal to that of OF and further takes the maximal value. 



Let X\ = I sm9 sin tp, X2 = Icosip and x 3 = ~ cos 9 sin ip with ^ 9 < 2n and 
^ ip < ir. It is not difficult to see that either of the following conditions give rise to 

W + \ = W-\- 

sin cos a — cos sin a = 0, (58) 

q 2 cos 2 9 - (1 - q) 2 cos 2 (9 - 2/3) = 2g - 1, (59) 

The condition (58) gives us a trivial result, namely, y + =y_=r B . So we focus our 
attention on condition (59). By lengthy but straightforward calculations, we find that 
under (59), the maximal value of |y+| or |y_| is given by 



2/max = y 1 - 4g(l - q) cos 2 a sin /3. 
Under this decomposition, we obtain the average entropy 

&EF — tl 



(60) 



(61) 



Combining (55), (60) and (61), we see that S m { n — Sef- Thus we attribute a geometric 
interpretation to the analytical expression for S'min. 



6. Two-qubit states with rank larger than two 

If r&nk(p AB ) > 2, the complement state p BC is a 2®r state with r = 3 or 4. Although 
ra.nk(p BC ) = 2, we cannot ensure that E(p BC ) is obtained on a two-element ensemble. 
In this situation, if we performed von Neumann measurement on A and acquire the 
MAE about qubit B, the MAE obtained in this way would not be the true EoF of 
p BC , but rather an upper bound for E(p BC ). 

Despite the fact that we have no definite results about E(p BC ), the upper and the 
lower bound are known. In [37], the /-tangle r, an entanglement measure proposed 



Quantum discord of two-qubit rank-two states 18 

by Rungta et al in [38], is calculated exactly for any bipartite states with rank two. 
For the state p BC under consideration, the following inequality holds [37]. 

E(P BC ) < h(± + yi-r(pBC)). (62) 

On the other hand, the exact concurrence of p BC is available from the approach 
presented in [27,28] (see also the Theorem 1). The following inequality is from [28]. 

E(P BC ) > h(± + yi-Con 2 (p*C)). (63) 

Recall the K-W relation S m in = E(p BC ). The inequality (62) and (63) give the upper 
and lower bound for the MAE respectively. In [39] , the authors provide the upper and 
lower bound for the quantum discord from this consideration. 

Below we discuss two issues. One will concern the upper bounds for the EoF of 
the rank-two state p BC , and the other the more general measurements performed on 
qubit A to achieve the MAE S'min. 

Suppose that we are restricted to make von Neumann measurements on qubit A. 
Let S min be the MAE of the postmeasurement ensemble for p B . It is an upper bound 

for E(p BC ), i.e., E(p BC ) < S min . Note that the right-hand side of the inequality 
(62), denoted by H(t) for short, is another upper bound for E(p BC ). It is claimed 
that the expressions on the left-hand and right-hand side of (62) usually differ only 
by about 10 . Let us compare these two upper bound. We select randomly X 

(2) 

states with rank three or four, calculate S min and H(t), and plot the difference 

—(2) (2) 

A = H(t) — S min in figure 4. In calculating S min , we adopt the approach provided 
in [40]. That is, the optimization only concerns the equi-entropy decomposition and 
quasi-eigendecomposition of p B [26] . It should be mentioned that this approach is not 
strictly correct [41-43]. 

—(2) 

In figure 4 we see that A Js 0, which means that S mm is a tighter upper bound 
for E(p BC ). For most of X states, A is about 10~ 15 . In fact the two upper bounds 
are equal to each other for these states. However there are some state for which the 
difference A can be about 10 -2 , much larger than 10 -4 . Further study shows that, 

(i) if S min is attained from the equi-entropy decomposition of p B , then A = 0; (ii) if 
S m i n is attained from the quasi-eigendecomposition, then the difference A is obviously 
larger than zero. 

For general two-qubit states, we calculate numerically the difference A for about 
10 4 randomly selected states and plot the results in figure 5. There are more states 
for which A > 0, meaning that the MAE is indeed a tighter upper bound for the EoF 
of the complement state. 

Now let us turn to considering the general POVM measurements on qubit 
A. In [28], the authors studied the entanglement entropy for axially symmetric 
channel. It has been found that in a small region of the parameter space the optimal 
decompositions are of length three. To put it in the case we are considering, it means 
that the EoF of some p B can be obtained on a 3-state ensemble. It immediately 
follows that the optimal ensemble for p B in the sense of the MAE (or quantum 
discord) is also a 3-state one, and that the optimal measurement on qubit A is a 
three-element POVM. In fact, while performing numerical computation, we find the 
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Figure 4. For about 10 4 randomly selected X states p AB with rank three or four, 

we calculate the MAE S min for p and the H(t) for the complement states p BC . 

— (2) 
The difference A = H(t) — S min is plotted in the figure. 




Figure 5. The numerical results of A for about 10 randomly selected 2-qubit 
states. 



„AB 



following example. 

/ 0.7 0.2795 \ 



0.15 

\ 0.2795 0.15 / 

It is a rank-three X state. If we perform von Neumann measurement on qubit A, the 
MAE can be calculated as 



(64) 



S { ^ m = 0.295127, 



(65) 

In this case, the optimal measurement operators are i(l ± a x ), which give rise to the 
equi-entropy decomposition of p B . 
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Now we perform a three-element POVM measurement on qubit A. The three 
measurement directions are given by 

fei = (0.929301, 0, -0.369322), 

k 2 = (-0.929301, 0, -0.369322), 

fe 3 = (0, 0, 1). 
The corresponding probability is 

Probi = Prob 2 = 0.365144, Prob 3 = 0.269712. 
Under such measurement the MAE is 

S^ n = 0.291942 < Sil- (66) 

(2) (3) 

Although the difference S min — S min w 0.003 is small, it indicates that the three- 
element POVM measurement is optimal. 

7. Conclusion 

In conclusion, we consider the K-W relation from the viewpoint of quantum channel. 
The quantum channel, isomorphic to the given state, determines the form of the 
quantum steering ellipsoid, which is a useful tool to discuss quantum discord. In 
the context of K-W relation, this channel is closely related to the concurrence of the 
complement state. 

Recently, several works are devoted to find more rigorous results of quantum 
discord for 2-qubit states [41-44]. Here we have shown that K-W relation, together 
with known results of EoF, provide us useful information for this problem. We point 
out that the quantum discord of any 2-qubit rank-two state can be evaluated exactly 
and strictly. In this case, the optimal measurement on one qubit is the von Neumann 
measurement, which will induce the equi-entropy decomposition for the reduced state 
of the other qubit. The situation is more complicated for the 2-qubit states with rank 
larger than two. Von Neumann measurement may not be the optimal choice. One has 
to take three-element POVM measurements into consideration. 
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